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Abstract 

In this paper we study the breakdown of normal hyperbolicity and its consequences for reaction 
dynamics; in particular, the dividing surface, the flux through the dividing surface (DS), and the 
gap time distribution. Our approach is to study these questions using simple, two degree-of-freedom 
Hamiltonian models where calculations for the different geometrical and dynamical quantities can 
be carried out exactly. For our examples, we show that resonances within the normally hyperbolic 
invariant manifold may, or may not, lead to a 'loss of normal hyperbolicity'. Moreover, we show 
that the onset of such resonances results in a change in topology of the dividing surface, but does 
not affect our ability to define a DS. The flux through the DS varies continuously with energy, 
even as the energy is varied in such a way that normal hyperbolicity is lost. For our examples the 
gap time distributions exhibit singularities at energies corresponding to the existence of homoclinic 
orbits in the DS, but these singularities are not associated with loss of normal hyperbolicity. 

PACS numbers: 05.45.-a, 45.10.Na, 82.20.Db, 82.20.-w 



I. INTRODUCTION 

Transition state theory has played, and continues to play, a fundamental role in how 
we think of chemical reactions. There are many excellent reviews of the subject; see, for 
example, refs [TTl2ll3Mlf5l The work of Wigner^ highlighted the notion of the transition state 
as a dividing surface (DS) in phase space, having the property that trajectories crossed this 
surface (only once) in their evolution from reactants to products. This phase space point 
of view was relatively undeveloped for many years until the 1970's when such a dynamical 
theory was fully realised for two degrees of freedom (DoF) in the work of Pechukas, Pollak 
and ChilcPSMMiilMil] j n their work they realised the construction of such a DS having 
the 'no recrossing' property, and a central role was played by a saddle-like (in terms of 
stability) periodic orbit (PO) which serves as the 'anchor' for the definition of the DS. This 
PO was called a periodic orbit dividing surface (PODS). (This choice of name was perhaps 
unfortunate as the PO itself is not the DS [in phase space] but rather the boundary of the 
DS.) The work of Pechukas, Pollak and Child just cited was the first example of the use of 
a normally hyperbolic invariant manifold (NHIM) in the context of transition state theory 
(although see also ref. HM . 

For more than two DoF a suitable generalisation of a saddle like (or 'hyperbolic') periodic 
orbit was required, and this was realised in the notion of a normally hyperbolic invariant 
manifold (NHIM). Explicit discussion of NHIMs in the context of reaction theory was in 
refs [TBTfTBI (see also refs [TTffTSl) . However, the usefulness of NHIMs for the definition of phase 
space dividing surfaces was fully realised only after the development of methods for com- 
puting NHIMS and their associated DS. In particular, NHIMs could be computed using the 
classical Poincare-Birkhoff normal form procedure, and its quantum analog^ 20 * 21 *. However, 
the Poincare-Birkhoff normal form procedure is local in nature since it involves a Taylor 
expansion about an appropriate saddle point. Practically, the expansion must be truncated 
at an appropriate order so that its ability to describe geometrical structures and dynamics 
at a desired accuracy must be assessed (and this is, typically, problem dependent). Very 
broadly speaking, the NHIM, and its associated local dynamics, are accurately described by 
the truncated normal form for energies 'close' to the energy of the saddle point. To date, 
our knowledge of NHIMs for specific problems relies heavily on the techniques used to com- 
pute them, in particular normal form theory. The local nature of these techniques makes 



NHIMs tools of limited utility in terms of our ability describe reaction dynamics at energies 
well above the reaction threshold. In particular, very little is known about how the NHIM 
behaves as we increase the energy above that of the saddle. We expect the normal form 
to 'break down' (in some not very precisely defined sense). However, the breakdown of the 
normal form and the breakdown of the NHIM are two separate issues. Questions concerning 
the existence and properties of the NHIM can in principle be considered independently of 
the methods used to compute the NHIM, when normal form approximations may no longer 
apply. 

In this paper we study the behaviour of NHIMs as the energy above the saddle point 
is increased, and the dependence on energy of associated quantities related to reaction 
dynamics. Our strategy will be to avoid the use of normal forms by considering simple 
models that allow us to explicitly compute the NHIM, the DS, and the geometrical properties 
associated with reaction dynamics. We can view our approach as a decoupling of normal 
form issues from the dynamics of reaction. We then ask the question: "what happens to a 
NHIM as parameters (the energy being a very important parameter) are varied"? In other 
words, we study 'bifurcation of NHIMs'. But, more specifically for reaction dynamics, we 
want to understand how such bifurcation of NHIMs affects the various quantities (described 
below) used to describe reaction dynamics. 

First, we need to describe some background regarding what we mean by the 'bifurcation 
of a NHIM' and then discuss our approach for studying these phenomena. Roughly speaking 
(for details, see Section [H]) a normally hyperbolic invariant manifold has the property that 
the linearized growth rates normal to the manifold dominate the linearized growth rates 
tangent to the manifold. The nature of the dynamics within the invariant manifold can 
be arbitrary, provided the growth rate conditions are satisfied. Hence, there are several 
possibilities to consider: 

• Bifurcation within the NHIM, but with the normal hyperbolicity conditions not af- 
fected. 

• Breakdown of the normal hyperbolicity conditions. 

• A combination of the above two phenomena - bifurcation within the NHIM and break- 
down of the normal hyperbolicity conditions. 



How bifurcation within the NHIM occurs and how the normal hyperbolicity conditions 



can break down will be discussed in detail when we consider explicit models in Section HI A 
Briefly, there are two ways in which such bifurcations can occur. One is when a parameter 
in the Hamiltonian is varied, and the other is when the energy itself is varied. (We will be 
considering only Hamiltonian systems.) Of course, energy is a parameter, but its importance 
in understanding the dynamics is such that we highlight its effects explicitly. 

Typically in bifurcation theory simple 'normal forms' are developed that embody the 
essential features of a particular bifurcation. We will follow a similar approach here. Our 
models will be two degree-of-freedom, separable, Hamiltonian systems where we have ex- 
plicit control and understanding of the dynamics in such a way that we can model the three 
possibilities noted above. Modelling the three possibilities for separable, two DoF Hamil- 
tonian systems is of course fairly straightforward, but that is the point. Our goal is to 
determine the effect that the three scenarios outlined above have on quantities that are used 
to describe reaction dynamics. In particular, we want to consider the effect on: 

• The dividing surface (DS) separating reactants and products as a function of energy 



(the DS is defined and explained in Section III B). 

The flux through the DS as a function of energy. 

The gap time distribution as a function of energy (these quantities are defined and 



explained in Section II B). 

Despite the fact that our examples are very simple (completely integrable 2 DoF Hamil- 
tonian systems), they allow for exact calculations and therefore provide benchmarks against 
which more complex examples can be compared. In particular, we note a previous study by 
Li and coworkers 22 which suggests that the breakdown of normal hyperbolicity of the NHIM 
is responsible for the non definability of a relevant DS at high energy above the reaction 
threshold. This conclusion is contrary to that reached in the study of our model examples 
(see below). A study of Allahem and B art sch 2 ^ asserts that the transition state loses, and 
then regains normal hyperbolicity as energy is varied. A study of Inarrea et aP^ concludes 
that even though bifurcation of the NHIM can occur, dividing surfaces can still be defined. 
These three studies are concerned with Hamilltonian models that are much more complex 
than ours (although still 2 DoF). We also note an interesting study by Yang 2 ^, although the 
models studied there are not Hamiltonian. 



The organisation of this paper is as follow. We review the concepts of NHIMs in Section 

we 



[H] where we discuss the definition of normal hyperbolicity in detail. In Section |II A 



discuss the general form for the two DoF models that we study. These models allow us 
to explicitly compute NHIMs and quantify their geometrical and stability properties. In 



Section |HB| we discuss concepts of phase space volumes, gap times, and reaction rates that 
are relevant for the present study. In Section HI A we describe the two models for which 



we will carry out explicit calculations. In Section III B we discuss dividing surfaces and 



compute the flux through the dividing surfaces for our two examples, while Section |IHC 
deals with the computation of the gap time distribution. In Section [III D we discuss the loss 
of normal hyperbolicity and its consequences for the DS, flux, and gap time distribution in 
our examples. In Section [TV] we present our conclusions. 



II. NORMALLY HYPERBOLIC INVARIANT MANIFOLDS AND MODELS FOR 
THEIR BIFURCATION 

The notion of a normally hyperbolic invariant manifold (or NHIM) is by now a standard 
concept and tool in dynamical systems theory that is playing an increasingly important role 
in a variety of applications. The theoretical framework developed over the course of many 
years, beginning in the early part of the 20th century and reached a mature form in the 
works of Fenichet 26 * 27 * 28 ! and Hirsch, Pugh, and ShutPH In ref. EHl a (relatively) elementary 
exposition of Fenichel's approach to NHIMs is given with some discussion of applications. 

We begin by defining the notion of a normally hyperbolic invariant manifold (NHIM) in 
a continuous time setting (i.e. the setting of autonomous ordinary differential equations, 
rather than maps). There are numerous definitions that appear throughout the literature. 
We take the definition from ref. [311 From an 'applied' point of view the definition may 
appear somewhat abstract, leading to difficulties understanding how the definition can be 
applied in concrete settings. However, this level of abstraction actually provides a great deal 
of flexibility for applications, as we shall see when we consider explicit models for illustrating 
different types of bifurcation phenomena associated with NHIMs and reaction dynamics. 

It is sufficient for our purpose to take phase space to be M 2n , which is a C°° differentiate 
manifold. For inner products and norms on this space we will take the usual Euclidean 
structures. Phase space is taken to be even dimensional since we will be considering canonical 
Hamiltonian systems, i.e. we consider equations of the following form: 

BH > < • - (2.1a) 

(2.1b) 

for some function H(qi, . . . , q n ,Pi, ■ ■ ■ ,Pn) = H(q,p). We require the Hamiltonian H(q,p) 
to be at least C k+1 since we will assume the flow that it generates is at least C k . In 
practice, this is not a restriction since the explicit Hamiltonians that we consider will be 
infinitely different iable. However, in the theoretical description of NHIMs, it is important 
to explicitly specify the degree of differentiability of all geometrical objects for reasons that 
will become apparent as we develop the discussion of NHIMs further. Note that in this 
work we will have no need to consider IR 2n as a symplectic vector space, or to investigate the 
associated consequences of symplecticity; this is likely to be a very fruitful avenue for future 
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research. The reason is the following. Roughly, in Hamiltonian dynamical systems stability 
properties are divided into 'elliptic' and 'hyperbolic', where bifurcation may occur when 
passing through the boundary between stability types. Results describing 'elliptic stability' 
(e.g. the Kolmogorov-Arnold-Moser (KAM) and Nekhoroshev theorems) rely heavily on 
the Hamiltonian structure for their proof. 'Hyperbolic stability' results, on the other hand, 
do not generally rely on the Hamiltonian structure, and that is true for the development of 
NHIMs. It would nevertheless be interesting to know what restrictions Hamiltonian structure 
places on the structure of NHIMs. For example, in ref . E21 it is proven that invariant tori in 
canonical, time independent, Hamiltonian systems cannot be normally hyperbolic. 

We suppose that $* : M 2n x R — > R 2n is a C k smooth flow defined on R 2n generated by 



(2.1 ). We will also need assume the flow to exist for all time, i.e. t G R. As phase space is 
not compact, we will need to explicitly discuss the issue of existence of the flow for all time 
in our examples. 

We now state the definition of a NHIM sufficient for our needs. 

Definition 1 (NHIM) An l dimensional submanifold, A of R 2n (£ < 2n) is said to be a 
normally hyperbolic invariant manifold for $* if A is invariant under $ 4 and there exists a 
splitting of the tangent bundle of TR 2n into sub-bundles: 

TR 2n = E u © E s © TA, (2.2) 

having the following properties: 

Invariance. The sub-bundles are invariant under D& , for all x G A, t G R, i.e., 

veE s x ^ D¥(x)(v) C E% t{x) for all x G A, teR, (2.3a) 

veE%^D&(x)(v)cE% tix) for all x E A, teR, (2.3b) 

v eTA x ^ D&(x){v) cTA^t^) for all x G A, teR (2.3c) 

Growth Rates. There exists a constant C > and rates < (5 < a such that for all x G A 
we have: 

veE 8 x =>\\ D¥(x){v) \\< Ce- at || v || for all t > 0, (2.4a) 

veE%=>\\ D¥(x)(v) \\< Ce at || v || for all t < 0, (2.4b) 

v gT^A^H D¥(x)(v) \\< Ce m || v || for all teR, (2.4c) 
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We now provide some additional background related to the concepts in this definition. 

Invariance. Invariance of A under $* means that for any x G A, $*(#) G A for all iel. 
We have not required A to be boundaryless. Making the assumption that A has no 
boundary is usual when considering invariance of A, since otherwise one would need to 
make some sort of assumption preventing trajectories from leaving A by crossing the 
boundary of A. We have not required A to have no boundary since that is not natural 
for the examples we will consider. Since we will be considering only autonomous 
Hamiltonian systems the fact that trajectories are prevented from leaving by crossing 
the boundary of A is ensured by constancy of the Hamiltonian function, i.e. energy 
conservation. 



Rates. Note the rates < (3 < a. These, along with (2.4a), (2.4b), and (2.4c) encapsulate 
the idea of normal hyperbolicity, in the sense that the rate of growth or decay of 
tangent vectors normal to A under the linearized dynamics dominate the rate of growth 
or decay of vectors tangent to A under the linearized dynamics. Note that vectors 
tangent to A can still exhibit exponential growth or decay, but the linearized rate 
must be smaller than that normal to A. 

Energy Conservation. The level set of the Hamiltonian, 

E E = {(q,p)eR 2n \H(q,p) = E}, (2.5) 

is a In — 1 dimensional surface (except at possible singular points of the surface) 
that is invariant under the flow generated by Hamilton's equations. If it is compact 
and boundaryless, then the flow exists for all time. Generically 33 the m dimensional 



manifold A will intersect (2.5 ) in a £—1 dimensional manifold, which we refer to as Ag. 
If A is normally hyperbolic, then A^ = A fl £# is also normally hyperbolic. Thus, for 
Hamiltonian systems we will be considering both iso-energetic and non iso-energetic 
NHIMs. 

Dimensionality. We will assume that: 

dim^ = dim£* = m, (2.6) 

and it follows from the conditions of Definition IT] that T X A, E x and E s x vary continu- 
ously with respect to x, and therefore the dimensions of T X A, E x and E x are constant 



with respect to x. We also have: 

&imT x k = L (2.7) 



It follows from (2.2) that: 



2n = 2m + L (2.8) 



Hence, I is even. 



The Need for Bundles: Linearization Characteristics. The 'bundle' concept and 
terminology, while perhaps unfamiliar, is necessary in our setting since we are con- 
sidering the linearization about an invariant manifold, which is generally filled with 
different orbits (as opposed to the more familiar situation of linearizing about a sin- 
gle orbit, such as an equilibrium point or periodic orbit). Consider the unstable 
bundle, for example. Invariance of this bundle under the linearized flow means that 
v G E% =$> D^ t (x)(v) C E^ t ( x ) for all x G A, t G M.. This relation provides an 
unstable growth direction all along the orbit through x, i.e. all along $ 4 (x), iel, and 
enables us to quantify the unstable growth along a given orbit. The (disjoint) union 
of all of the unstable subspaces over all points in A gives the unstable bundle over A. 
Similar consideration hold for the stable bundle over A and the tangent bundle of A. 
More details can be found in the general references for NHIMs given earlier. 

a. Existence of Stable and Unstable Manifolds. As we have previously mentioned, the 
definition of NHIM given in Definition [T] specifies conditions on the linearized dynamics about 
A. These provide sufficient conditions for proving results about the nonlinear dynamics. In 
particular, the sets of points having the following properties: 

W S {A) = {y G R 2n | d ($%, M) < C y e~ at for all t > 0} (2.9a) 

W U (A) = {y G R 2n | d ($*(y, M) < C y e at for all t < 0} (2.9b) 

can be shown to be C s ~ l manifolds (for some constant C y > 0), where s < min < k, § >, where, 
we recall, k is the degree of differentiability of the flow. Moreover, under the conditions 
of Definition [TJ A can be shown to be a C s invariant manifold. However, we emphasise 
again that the conditions of Definition [I] provide sufficient conditions (in terms of linearized 
dynamics) for the smoothness properties of the manifolds. For example, in the models that 
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we consider in Section HI A A is explicitly known and is infinitely differentiable, regardless 



of the value of the ratio |. Similarly, the flow will be infinitely differentiable in the models 
that we consider in Section IHI A[ 

b. Persistence Under Perturbation. One of the most important properties of NHIMs 
is that they persist under C l perturbations. It is important to clearly define the term 
'perturbation'. For example, we are considering time independent Hamiltonian systems. 
Time-dependent erturbations would not satisfy the hypotheses of the perturbation theorem, 
unless the time-dependence could be recast in a way that the system became time indepen- 
dent, and the resulting system satisfies the hypotheses of the perturbation theorem. Time 
dependence that is periodic or quasi periodic can be treated in this way, while more general 
time dependence results in a breakdown of compactness, leading to problems satisfying the 
hypotheses of the theorem. 

c. Sufficient Conditions and 'Breakdown of Normal Hyperbolicity ' Definition [I] is a 
definition of a NHIM. It cannot be concluded that if the conditions of the definition are not 
satisfied, such as might occur for certain parameter values in a parameterized system, that 
the NHIM undergoes a 'bifurcation'. Loss of hyperbolicity is a necessary, but not sufficient 
condition for bifurcation to occur, but the nature of a bifurcation depends on nonlinearity, 
which is outside the standard requirements of normal hyperbolicity as stated in Definition 

m 

A. The General Form of the Models Under Consideration 

We now consider the general form of a model that encompasses all of the specific examples 
that we will study. We consider a two DoF Hamiltonian system of the following form: 

H(qi,pi,q2,p 2 ) = H 1 (qi,pi) + H 2 (q2,p 2 ), (qi,Pi, ^2,^2) e K 4 . (2.10) 
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with associated Hamiltonian vector field: 



?i = -~-^(qi,Pi), (2.11a) 

dpi 

f) IT 

Pi = ~(quPi), (2.Hb) 

dqi 

& = i^(g2,P2), (2.11c) 

dp 2 

f) IT 

P2 = - 7 r 1 (q2,P2). (2.11d) 

dq 2 



Note the particularly simple form of eqs (2.11 ). It has the form of two separable one degree 



of-freedom Hamiltonian systems (hence it is completely integrable). The flow generated by 



(2.11) has the general form: 

$*(<?io,Pio,920,P2o) = (qi(t,qio,Pio),Pi(t,qio,Pw),q2(t,q20,Pm),P2(t,q20,P2o))- (2-12) 
The equation for the (typically) three-dimensional energy surface is given by: 

^e = {(qi,Pi,q2,P2) e M A \H(q 1 ,p 1 ,q 2 ,p 2 ) = H x (q\,p\) + H 2 (q 2 ,p 2 ) = E) . (2.13) 



We make the following assumption on the q 2 — p 2 component of (2.11). 



Assumption 1 At q 2 = p 2 = the system: 

PjTT 

g 2 = ^(g 2 ,p 2 ), (2.14a) 

op 2 

P2 = -^(q2,P2). (2.14b) 

oq 2 

has a hyperbolic equilibrium point. The (positive) eigenvalue of the matrix associated with 



the linearization of (2.14) about this equilibrium is a > 0. 

Clearly, the set 

A = {(qi-Pi,q2,P2)\q2 = P2 = 0} , (2.15) 

is a two dimensional invariant manifold in the four dimensional phase space. The intersection 
of the 3D energy surface with this 2D invariant manifold is given by: 

A B = AnE E = {(qi,Pi,q2,P2) e M 4 | H( gi , Pl ,0,0) = H 1 (q 1 ,p 1 ) + H 2 (0,0) = E) , (2.16) 

which is (typically) a one dimensional level set of H\(q\,pi)\ an isoenergetic invariant man- 
ifold. Now we will show that, under Assumption [TJ A and A^ are both normally hyperbolic 
invariant manifolds. We therefore need to show that Definition [Q holds for A. 
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We begin by computing the linearization of the flow of (2.12) about an arbitrary point on 



A, which is denoted by x = (gi,Pi, 0, 0). First, we transform -f^fe, P2) to a set of coordinates 
that facilitates the computations. It follows from Assumption [T] that there exists a linear, 
symplectic transformation of coordinates, (<?2,£>2) — > (^2,^2), such that the quadratic part of 
H2 is "diagonal" in these coordinates, i.e. 



H 2 (q2,P2) = #2(0,0) + ap 2 q2 + H (q 2 ,P2 



(2.17) 



where H 3 (q 2 ,p2) is 0(3) (note: there are no linear terms in H 2 {q2,P2) since (^2,^2) = 
(^2,^2) = (0, 0) is an equilibrium point). In these coordinates we can write the Hamiltonian 



(2.10) in the form: 

H(qi,pi,q~2,p2) = H 1 (q ll p 1 ) + H 2 (q2,P2), (qi,Pi,q~2,p2) e M 4 . 



(2.18) 



In these coordinates, the linearisation of the flow has the the following block diagonal form: 



D¥(x) 



( 



\ 



A 



02x2 
,at 



\ 



J2x2 



e 



-at 



(2.19) 



/ 



where 



.4 



9 2 Hi („ „ \ d 2 H! 



(2.20) 

O2X2 denotes the 2x2 matrix of zeros, and we assume that we have transformed the ^2 ~P2 
coordinates of the lower right hand 2x2 block so that the flow assumes the diagonal form, 
which follows from Assumption [TJ 

A general tangent vector in E%, where we use the shorthand notation x = (51, pi, 0, 0), 
has the form 

/o\ 





1 



(2.21) 



Then we have: 



D&(x)(v) 



( \ 


e at 

\ } 



(2.22) 
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which is a tangent vector in E%,, where x' = $*(x). Since this argument holds for any x E A 



and any v E E™ it follows that the bundle E u is invariant (i.e. (2.3b) holds). Moreover 



it follows from (2.22) that the growth rate condition (2.4b) holds for all vectors in E u . A 



similar argument for the invariance of E s and the growth rate of vectors in E s follows. 

Now we turn our attention to the tangent bundle of A. A general vector in T X A has the 
form: 



b 





(2.23) 



Then using (2.19) we have: 



D¥(x)(v) 



fa'\ 
V 


V 0/ 



(2.24) 



which is in T X >A, where x' = $*(x). Here, a' and b' are functions of (qi,Pi), but the explicit 
functional form is not important for this argument. Since this argument holds for any x E A, 
v E T x /A, it follows that TA is invariant. What we can not verify at this point is the growth 



rate condition (2.4c). This will require explicit conditions on the flow on A. The flow on 



A can be very general, but normal hyperbolicity will require that (2.4c) is satisfied by this 



flow. This condition will be considered in the specific examples that we analyse in Section 

mm 

B. Phase space volumes, gap times, and reaction rates 

In this section we briefly review the concepts from classical reaction rate theory that are 
relevant to the present study. This section is adapted from the paper of Collins et aP^ where 
more background and details can be found. 

Points in the 4-dimensional system phase space M. = M 4 are denoted z = (pi,P2,Qi,Q2) = 
(PjQ) *= -M- The system Hamiltonian is denoted by H(z), and the 3 dimensional energy 
surface at energy E, H(z) = E, is denoted £# C M.. The corresponding microcanonical 
phase space density is S(E — H(z)), and the associated density of states for the complete 
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energy surface at energy E is 

p(E)= f dz5(E-H(z)). (2.25) 

Jm 

The disjoint regions of phase space ('reactant', q 2 < 0, and 'product', q 2 > 0) separated 
by the phase space dividing surface DS(-E) are denoted A4±; the region of phase space 
corresponding to q 2 > will be denoted by Ai+, and that corresponding q 2 < will be 
denoted by M... 

The microcanonical density of states for points in region Ai + is 

p + (E)= f dzS(E-H{z)) (2.26) 

Jm+ 

with a corresponding expression for the density of states p.{E) in Ai_. Since the flow is ev- 
erywhere transverse to DS±(2£), those phase points in the region .M+ that lie on trajectories 
that cross DS±(-B) can be specified uniquely by coordinates (p, q, ip), where (p, q) G DS+(i£) 
is a point on DS + (.E), specified by 2 coordinates (p,q), and ip is a time variable. The 
point z(p,q,ip) is reached by propagating the initial condition (p,q) G DS + (£ I ) forward 
for time ^l 35 l 36 l 37 l. As all initial conditions on DS + (.E) (apart from a set of trajectories of 
measure zero lying on stable manifolds) will leave the region A4 + in finite time by crossing 
DS.(-E'), for each (p,q) G DS + (£'), we can define the gap time s = s(p,q), which is the 
time it takes for the trajectory to traverse the region M. + before entering the region A4-. 
That is, z(p,q,if) = s(p,q)) G DS.(.E'). For the phase point z(p,q,ip), we therefore have 
< ij} < s(p,q). 

The coordinate transformation z — > (E,ip,p,q) is canonica F 5 * 38 * 39 * 4 ^ , so that the phase 
space volume element is 

d 4 2 = dEd^da (2.27) 

with der = dpdq an element of 2 dimensional area on DS(E). 

The magnitude 4>(E) of the flux through dividing surface DS+(.E7) at energy E ('direc- 
tional flux') is given by 



cj>{E) 



da 

T>$+(E) 



(2.28) 



where the element of area da is precisely the restriction to DS(E) of the appropriate flux 
2-form w corresponding to the Hamiltonian vector field associated with ff(^)[Mi i | 42 | 43 | rj^g 
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reactant phase space volume occupied by points initiated on the dividing surface with ener- 
gies between E and E + dE is therefore^ 11 39 1 40 1 44 1 45 1 46 1 47 1 



dE da dip =dE da s (2.29a) 

iDS+(E) JO JdS+(E) 



dE <f){E)s (2.29b) 



where the mean gap time s is defined as 

0(£) 



s = -jj^ I da s (2.30) 

DS+(£) 



and is a function of energy E. 



1. Gap time and reactant lifetime distributions 

The gap time distribution, V(s; E) is of central interest in unimolecular kinetics^^S the 
probability that a phase point on DS + (£') at energy E has a gap time between s and s + ds 
is equal to V(s; E)ds. An important idealized gap distribution is the random, exponential 
distribution 

V(s; E) = k(E) e- k{E)s (2.31) 

characterized by a single decay constant k (where k depends on energy E), with correspond- 
ing mean gap time s — k^ 1 . An exponential distribution of gap times is usually taken to be 
a necessary condition for 'statistical' behavior in unimolecular reaction d 35 * 48 * 49 * 50 -*^. 

The lifetime (time to cross the dividing surface DS.(-E')) of phase point z(p,q,ip) is t — 
s {Pi Q)~' l l ) i an d the corresponding (normalized) reactant lifetime distribution function P(t; E) 
at enerev E j j 351481491501521531541 

P(t ; E) = -— Prob(t > t'\ E) (2.32a) 



dt' 

i p+oo 



t'=t 



s .It 



dsV{s]E) (2.32b) 



where the fraction of interesting (reactive) phase points having lifetimes between t and t + dt 
is P(t; E)dt. It is often useful to work with the unnormalized lifetime distribution F, where 
F(t;E) = s¥(t;E). 
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Equation (2.32a) gives the general relation between the lifetime distribution and the frac- 



tion of trajectories having lifetimes greater than a certain value for arbitrary ensemble d 52 * 53 * 54 '. 



Note that an exponential gap time distribution (2.31) implies that the reactant lifetime dis- 
tribution P(t; E) is also exponentia l 35 ' 48 ' 49 * 52 * 53 * 54 !; both gap and lifetime distributions for 
realistic molecular potentials have been of great interest since the earliest days of trajectory 
simulations of unimolecular decay, and many examples of non-exponential lifetime distribu- 
tions have been found 36 * 52 * 53 * 54 * 55 * 56 * 57 ' 58 ' 59 '. 
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III. MODEL DYNAMICS: DIVIDING SURFACES, FLUX, BREAKDOWN OF 
NORMAL HYPERBOLICITY, AND GAP TIME DISTRIBUTIONS 

In this section we describe our two model systems. Each model has NHIMs, and we 
describe the normal hyperbolicity properties of the NHIMs as a function of energy and 
other parameters. We then compute the flux through the associated DS as well as the gap 
time distribution, and we discuss the behaviour of these quantities in relation to the property 
of normal hyperbolicity. 

A. The Models 

The two models to be studied are particular examples of the general class of models 
described in Section III Al 

1. Example 1. Uncoupled Simple Pendulum and Symmetric Double-Well 

Our first example will consist of an integrable system whose Hamiltonian is the sum of 
the Hamiltonian of a simple pendulum and the Hamiltonian of a symmetric double-well: 

2 2 i 

H = -± -a 2 1 cosq 1 + -j - -^-q\ + -q 2 = H x (q x , p x ) + H 2 (q 2 ,p 2 ), (3.1) 

with associated equations of motion: 

dH 

q x = — = p u (3.2a) 

dpi 

dH 2 . , . 

Pi = -7^— = -«i singi, (3.2b) 

dq x 

dH 

<?2 = -w- = P2, (3.2c) 

dp 2 

PjTT 

P2 = ~a— = "2^2 - <?2 ( 3 - 2d ) 

oq 2 
The phase space structure for these two, uncoupled, one degree-of-freedom Hamiltonian 

systems is illustrated in Fig. [11 Since the trajectories of (3.2) can be solved analytically, for 



any initial condition, it can be seen explicitly that the flow generated by (3.2) exists for all 
time. 



Following (2.15), the two dimensional non-isoenergetic surface: 

A = {(gi,Pi,<Z2,p 2 ) e M 4 | q 2 = V2 = 0} , (3.3) 
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is an invariant manifold for (3.6). For a 2 > 0, tangent vectors normal to A experience 
exponential growth and decay under the linearised dynamics. However, we cannot claim 
that it is a NHIM according to definition [TJ Whether or not this is in fact the case depends 



on the dynamics on A. Following (2.16), the intersection of the three dimensional energy 



surface (3.1) with (3.3) is given by: 



As = Uqi,pi, q 2 , p 2 ) <ER 4 \ -^ - a\ cos q 2 = E, q 2 = p 2 = j- . (3.4) 

Hence, the dynamics on A is described by the dynamics of a (symmetric) two well potential, 
and an orbit of the pendulum with energy E defines A#. We note that: 

• <7i = £>i = is an elliptic equilibrium point of the pendulum with energy E = —a\. 

• Surrounding the elliptic fixed point is a family of periodic orbits ('librations') whose 
energies increase monotonically from E = —ot\ to E — ot\. 

• q\ = it, pi = is a saddle point with energy E = a\ that is connected by a pair of 
homoclinic orbits. 

• Outside the pair of homoclinic orbits are two families of periodic orbits ('rotations') 
whose actions increase monotonically with energy, from E = a\. 

The phase space structure associated with the simple pendulum on A is shown in Fig. [TJ 
We note that homoclinic orbits have zero frequency. Hence, for one DoF Hamiltonian 
systems, they are examples of the phenomenon of resonance. The presence of homoclinic 
orbits therefore allows us to examine the role of resonance in the breakdown of normal 
hyperbolicity explicitly. 

In our model we take 'reaction' to correspond to a change of sign of the coordinate q 2 . 
Hence, the energy of reacting trajectories in the q 2 = p 2 coordinates must be greater than 
that of the saddle point at q 2 = p 2 = 0, i.e., the energy must be greater than zero. Hence, 
the total energy must be such that the energy in the q 2 —p 2 subsystem is greater than zero. 
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2. Example 2. Two Uncoupled Symmetric Double-Wells 



The second example consists of two uncoupled symmetric double-well Hamiltonians: 



H 



P\ Qfl 2 , 1 4 , V\ "2 2 . 1 4 

#l(?l,Pl) + H 2 (q 2 ,p 2 ), 



with associated Hamiltonian vector field: 



dH 



Pi 



dpi 
dH 



dqi 
dH 

q.2 = -7T- = Pi 



«1?1 -ffi! 



P2 



<9p 2 
5g 2 



«2g2 -g 2 - 



(3.5a) 
(3.5b) 

(3.6a) 
(3.6b) 
(3.6c) 
(3.6d) 



Again, the trajectories of (3.6) can be solved analytically for any initial condition, and the 



flow generated by (3.6) exists for all time. 



The discussion of the phase space structure is similar to that for Example 1 above. 



Following (2.15), the two dimensional non-isoenergetic surface: 

A= {{Qi,Pi,q2,P2) eR A \q 2 = P2 = o} 



(3.7) 



is an invariant manifold for ( |3.6[ ). For a 2 > 0, tangent vectors normal to A experience 
exponential growth and decay under the linearised dynamics. Again, whether or not A is 
a NHIM according to definition IT] depends on the dynamics on A. Following (2.16), the 



intersection of the three dimensional energy surface (3.5) with (3.7) is given by: 



A/; = { (qi,Pi, q 2 ,P2) G ^ I y - ^q\ + \ q\ = E, q 2 = p 2 = }. 



(3.8) 



Hence, the dynamics on A is described by the dynamics of a symmetric two well potential, 
and an orbit of the symmetric two well potential with energy E defines A#. We note that: 

• qi = Pi = is & saddle point of the symmetric two well potential with energy E = 0. 



qi = ±Joii, pi = are elliptic equilibrium points with energy 



■H 



Surrounding each elliptic equilibrium point is a family of periodic orbits whose energies 
increase monotonically away from the equilibrium point. At E = the periodic orbits 
merge into a pair of homoclinic orbits that connect the saddle point at the origin. 
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• As the energy increases from zero there is a family of (symmetric) periodic orbits that 
surround both potential wells whose actions increase monotonically with energy. 

The phase space structure associated with the symmetric double well on A is shown in 
Fig. U 

As in the previous model, 'reaction' is associated with a change of sign of the q<i coordinate. 
Hence, the energy of reacting trajectories in the q<i — P2 coordinates must be greater than 
the energy of the saddle point at q2 = P2 = 0, i.e. the energy must be greater than zero. 
The total energy must therefore be such that the energy in the q2 — P2 component of the 
model is greater than zero. 

B. Dividing Surfaces and Flux 

For each of our model Hamiltonians we can construct a dividing surface in phase space 
having the (local) no-recrossing property. The codimension one non-isoenergetic surface 
defined by q 2 = (locally) divides the phase space into two regions: one associated with 
the regions defined by ^2 > and the other associated with the region defined by ^2 < 0. 
The dividing surface restricted to a three dimensional fixed energy surface H (qi,pi, q 2 , P2) = 
Hi (qi,Pi) + H 2 (q2,P2) = E is given by: 

DS(E) = Uq 1 ,p 1 ,q,p 2 ) \q 2 = 0,H = H 1 (q 1 ,p 1 ) + ^ = e\. (3.9) 

This dividing surface has two halves: 

DS+(E) = l{qi, P i,q 2 ,p 2 ) \q2 = 0,H = H 1 {q 1 ,p 1 ) + ^ = E, p 2 > o} (3.10a) 

and 

DS-(E) = Uqi,Pi,q2,P2) \q 2 = 0,H = H 1 (q 1 ,p 1 ) + V -±p\ = E, p 2 < o} . (3.10b) 

These two halves meet at an invariant manifold: 

Ae = {{qi,Pi, q 2 ,P2) \q 2 = 0,H = H 1 {q 1 ,p 1 ) = E,p 2 = 0}. (3.11) 

The nature of this invariant manifold (especially the normal hyperbolicity property) depends 
on E, the dynamics on A#, and a comparison of the linearised growth rates normal and 
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tangent to E, as we will see for our examples. In particular, we will examine the inter- 
relation between these characteristics. 

We now describe conditions under which DS + (E) and DS^(E) are surfaces having the no 
(local) re-crossing property (and we will verify these conditions in our particular examples). 
These surfaces are defined by q2 = 0. Therefore points on these surfaces leave if: 

q2=P2^0, (3.12) 



and therefore it follows immediately from their definition that (3.10a) and (3.10b) have the 
(local) 'no-recrossing' property. 

We denote the directional flux across these hemispheres by 4>+(E) and <p-(E), respectively, 
and note that <j>+{E) + <j).(E) = 0. The magnitude of the flux is \<j>+(E)\ = \<j).(E)\ = <j>(E). 
The magnitude of the flux and related quantities are central to the theory of reaction rates, 
as we described in Section III BL 

1. Example 1. Simple pendulum plus Symmetric Double-Well 



i. Dividing surface. In Figure k^ we plot the dividing surface (i.e. (3.9)) for three 



different energies: E = 0.5 (energy below the separatrices of the pendulum), E = a\ = 0.8 2 
(energy equal to the energy of the separatrices of the pendulum) , and E = 1 (energy larger 
than the energy of the separatrices of the pendulum). We see that the DS undergoes a 
bifurcation as the energy passes through the energy of the separatrix. A central question 
is whether or not this bifurcation has any effect on quantities that are important for 
quantifying reaction dynamics, such as flux and gap times. 

b. Flux across the DS. By Stokes theorem, the directional flux across half the DS, at 
energy E, is the area enclosed by the invariant manifold A# on the (qi, p{) plane. This area 
is just the integral A = § pi dq\ and is related to the action variable of the simple pendulum 
by J = y~ § pi dqi so that we have the relation A = 2nJ. For the simple pendulum there are 
three different cases to consider for the calculation of the action integral: the librations, the 
separatrices, and the rotations. The computation of these integrals for the different cases is 
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reviewed in the appendix [A] Here we just note the results: 

J t (E) = — fe(M) + k?Ki(k) - Kxih)] (3.13a) 

7T 

J S {E) = ^i (3.13b) 

TV 

J r (E) = P^S 1 (k r ) (3.13c) 

K r TT 

where the subscripts /, s and r denote libration, separatrix and rotation, respectively. /Ci 
and Si denote the complete elliptic integrals of the first and second kind, respectively and 
ki and k r are the moduli of the complete elliptic integrals, and are functions of the energy 
E: 

1 / CM 

(3.14a) 
(3.14b) 




E + a( 



With these expressions in hand we obtain expressions for the directional flux <f)(E) across 
half the DS: 

cj)i{E) = 16a! [Exih) + kflCxik) - K x {ki)\ (3.15a) 

(f> a (E) = 16«i (3.15b) 

<f) r {E) = -^E 1 (k r ) (3.15c) 

For the separatrix and the rotation cases an extra factor 2 appears because one must count 
the area enclosed by the curves having positive and negative momentum. The graph of the 
flux (j)(E) is shown in Figure El 

2. Example 2. Uncoupled Symmetric Double-Wells 



Dividing surface. In Figure [5J we plot the dividing surface (i.e., the surface (3.9)) 



a 

for three different energies: E = —0.5 (energy below the separatrices of the two well 
potential), E = (the energy of the separatrices of the two well potential), and E = 1 
(energy larger than the energy of the separatrices of the two well potential). As for example 
1, the DS undergoes a bifurcation as the energy passes through the energy of the separatrix. 
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b. Flux across the DS. As for Example 1, the directional flux across half the DS is 
the area enclosed by the invariant manifold A# on the (gi, p\) plane. This area is the 
integral A = § p\ dq\ and is related to the action variable of the symmetric double- well by 
A = 2tt J. For the symmetric double- well there are three different cases for the calculation 
of the action integral depending on the energy being less, greater or equal to the energy of 
the saddle point. The integrals for the different cases are (cf. appendix [B| : 



ME) = -|^ [(2 - kfi&fa) + (2*? - 2)/C 1 (* 1 )] (3.16a) 



3/2 

HE) = ^- (3.16b) 

3n 



ME) = ^ [(2kl - 1)^ 3 ) + (1 - kDK^h)} (3.16c) 

where the subscript 1, 2 and 3 refers to the three different regimes of energy, that is, energy 
below, and above the energy of the saddle point, respectively. In theses expressions a and b 
are roots of the equation pi = (cf. Appendix [BJ) and 5 = b 2 — a 2 . All these 3 quantities 
are functions of the energy. As for the case of the pendulum, the moduli k\ and k% of the 
elliptic integrals are functions of the energy: 

(3.17a) 

(3.17b) 

(3.17c) 

m ■ (3 ' 17d) 

The directional flux <p(E) across half of the DS is then given by: 

ME) = ^^ 2 [(2 - fc?)£i(*i) + (2*? - 2)£i(* 1 )] (3.18a) 

o 3/2 

ME) = ^- (3.18b) 





/ 


a(E) = 


= J ax - y a\ + AE 




1 


b(E) = 


= Jan + \Ja\ + AE 


h(E) -- 


/ 8(E) 

' V b 2 (E) 


h.JE) = 


lb 2 (E) 



>(E) = 4^ P k l - l ) £ ^) + 0-- k l)^^)\ ■ ( 3 - 18c ) 



3«3 



A plot of the flux (j)(E) as a function of energy E is shown in Figure [6 
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C. Gap Time Distributions 

In this subsection we compute the gap time distribution of our two examples. Before 
proceeding to the calculation of the gap time distribution, we need to determine the gap 
time for reactant trajectories. In our two examples the determination of the gap times of 
reacting orbits is straightforward since the gap time is just half of the period of the periodic 
orbits in the double- well degree-of-freedom (q2~P2) for which the energy is greater than the 
energy of the saddle point (q 2 = 0, p 2 = 0). The period of these periodic orbits as a function 
of the energy is given by: 

T(E 2 ) = -^L=K{h{E 2 % (3.19) 

where E 2 stands for energy in the independent second degree-of-freedom, and the total 
energy of the two, uncoupled one degree-of-freedom systems is denoted by E? = Ex + E 2 . 
The gap time as a function of the energy E 2 is then given by: 

s{E 2 ) = -J^L)C(k 3 (E 2 )), (3.20) 

For clarity here we will use the notations / for the action associated with the first degree- 
of-freedom (pendulum for example one and double-well for example 2) and J for the action 
associated with the second degree-of-freedom. Because of the direct relation between the 
energies of the different degrees-of-freedom and the actions associated with these degrees- 
of-freedom , at constant energy we have: 

J = J(E T ,I), (3.21) 

so that we can formally write the gap time as a function of the total energy and the action I: 
s(E 2 ) = s(Et, I). Again at constant Et, the number of trajectories having s < s = s(Et, I) 
is proportional to the area enclosed by the curve I — I in the (qi-pi) plane, that is to say, 
it equals 2ttI. The number of trajectories n for which the gap time s has the property that 

s < s < s + As is: 

dl 

n = 2n [(I + A/) - J] = 2tt— As. (3.22) 
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The gap time distribution is now given by: 

71 

P(s; E T ) = — (3.23a) 

dl , 

= 2tt— 3.23b 

OS 

= -£?-• (3-23c) 

ai\E T 

Figure [7] shows the gap time distribution for both examples. 

Figure [7] shows that the gap time distributions for both examples possess a singularity. 
This singularity arises from the vanishing of the derivative of the action / with respect to the 
gap time s, which is equivalent to the singularity of the derivative of this action with respect 
to the energy E\. This singularity appears because the action variable I is a piecewise defined 
function with respect to the energy E\. Whereas the action J is a continuous function with 
respect to the energy E\, its derivative (the inverse of the associated frequency) diverges at 
the separatrix (homoclinic orbits), which is exactly where bifurcation occurs. Interpreted 
geometrically, we see that just below the bifurcation the rate of growth of the area enclosed by 
the periodic orbits in the (qi-pi) plane increases to infinity as we get closer to the homoclinic 
orbit. In the same way, just after the bifurcation this growth rate decreases from infinity. 
The appearance of homoclinic orbits is manifested in the flux <J)(Et), which is directly related 
to the action J, through the appearance of an inflexion point at the energy of the bifurcation. 

D. Loss of normal hyperbolicity and its consequences on the dynamics of reaction. 

In this section we will turn our attention to the question of normal hyperbolicity of 
the isoenergetic invariant manifold, which forms the boundary of the isoenergetic dividing 
surface. For both examples A^ is either a periodic orbit or the union of homoclinic orbits 
and the saddle points that they connect. The two dimensional non-isoenergetic invariant 
manifold is defined by qi = P2 = 0, which is a hyperbolic saddle point with associated growth 
rate «2- For «2 > tangent vectors normal to A experience exponential growth and decay 
under the linearised dynamics. Hence one way that normal hyperbolicity could be lost is 
for «2 to go from positive to negative. However, this possibility is not interesting from the 
point of view of reaction dynamics since this would destroy the reactive (double well) nature 
of the problem. Hence, we will require a 2 > 0. 
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The more interesting situation is the violation of the growth condition (2.4c) in the 
definition of a NHIM (Definition flj. When A^ is a periodic orbit, this growth rate is zero, 
and A E is a NHIM. When A^ is a union of homoclinic orbits and the saddle points that they 
connect, whether or not A E satisfies the growth rate conditions of Definition [I] depends on 
the nature of a,\ and 02- In particular, for Example 1 ±ai are the eigenvalues associated 
with the saddle point on the q± — p\ plane and for Example 2 ± v /aT|" are the eigenvalues 
associated with the saddle point in the q± = pi plane. 

We therefore have the following situations: 

Example 1: For a\ > y / «2 the growth rate conditions of Definition flj are not satisfied, and 
A E is not a NHIM. 

Example 2: For <y.\ > a<i the growth rate conditions of Definition [I] are not satisfied, and 
A E is not a NHIM. 

We now examine the consequences of this 'loss of normal hyperbolicity' in more detail. 
1. Consequences of loss of normal hyperbolicity 



We imagine the parameters a± > a/o2 fixed for example 1 and a.\ > ct2 for example 2, 
and we vary the energy in such a way that we pass through the energy value corresponding 
to the homoclinic orbits and saddle points. With a.\ and 0:2 fixed as above, at this energy 
value the conditions of Definition Q] do not hold. This situation could be referred to as 'loss 
of normal hyperbolicity'. We examine the implication of this loss of normal hyperbolicity 
for the quantities that we have computed for our two examples. 

The Dividing Surface and the (Directional) Flux Through the Dividing Surface. 

The dividing surface, as a function of energy, is shown in Figure [3] for example 1 
and in Figure [5] for example 2. While in both examples the geometry of the surface 
undergoes a qualitative change is we pass through the bifurcation, we are still able to 
define a dividing surface having the no recrossing property as we pass through the 
bifurcation, i.e., even as the NHIM experiences a loss of normal hyperbolicity (as just 
defined) . 
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The directional flux, as a function of energy, is shown in Figure [4] for example 1 and 
in Figure [6] for example 2. For both examples we see that the flux varies continuously 
as a function of the energy, even as we pass through the bifurcation, i.e., even as the 
NHIM experiences a loss of normal hyperbolicity. 

The Gap Time Distribution. From Figure [7] we have already noted note that the gap 
time distributions for both examples possess a singularity. This singularity is associ- 
ated with the existence of a homoclinic orbit in the DS, and is not related to the loss 
of normal hyperbolicity of the NHIM. 
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IV. CONCLUSIONS AND OUTLOOK 

In this paper we have studied the breakdown of normal hyperbolicity and its consequences 
for quantities related to reaction dynamics; in particular, the dividing surface, the flux 
through the dividing surface, and the gap time distribution. 

Our approach is to study these questions using simple, two degree-of-freedom Hamiltonian 
models for which calculations for the different geometrical and dynamical quantities can 
be carried out exactly. For our examples, we showed that resonances (homoclinic orbits) 
within the normally hyperbolic invariant manifold may, or may not, lead to 'loss of normal 
hyperbolicity'. Moreover, we showed that for our examples the onset of such resonances 
results in a change in topology of the dividing surface, but it does not affect our ability 
to define a dividing surface (DS), and that the flux through the DS varies continuously 
with energy, even as the energy is varied in such a way that normal hyperbolicity is lost. 
For both our examples we have shown that the gap time distribution exhibits a signature 
singularity at energies corresponding to emergence of a homoclinic orbit in the DS, but these 
singularities are not associated with loss of normal hyperbolicity. 
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Appendix A: Action variables for the simple pendulum 

In this appendix we review the compution of action variables for the integrable simple 
pendulum system. The Hamiltonian of the simple pendulum is: 

H = — — a\cosq. (Al) 

In order to determine the action variable for this system we have to consider three different 
cases: libration (E < of), separatrix (E = af), and rotation (E > a\). 

1. Libration 

The action for this case is given by: 

Ji = ^jpdq. (A2) 

For this case the system doesn't have enough energy to cover the full range of the angle q 
and there should be two turning points where the momentum vanishes, p = 0. Let go be the 
positive value of q at the turning point which depends on the energy: 

go = cos -1 ( 2 ) • ( A3 ) 



The momentum p can be expressed as a function of go and g: 



p = \2a\(cosq — cosg ). (A4) 



Substituting this expression into the eq. (A2) we obtain the integral: 



4«i f qo / . 2 /g \ . 2 (q 



7r J V " \2/ V2 



J t = — Jsin 2 Hr - sin [7i)dq. (A5) 



After an appropriate change of variable we can put the former integral into a form which 
makes it resemble an elliptic integral: 

T 8a x kf r /2 cos 2 ,, lK . 



f Jo a/ 1 - k f sin2 v 
where the modulus ki of the elliptic integral is given by hi = sin (y). This integral can be 
expressed in terms of complete elliptic integrals as follows: 

J t = ** [S.ik) + k^k) - /d(fc,)] , (A7) 

7T L J 

where /Ci and E\ stand for the complete elliptic integrals of the first and second kinds 
respectively^. 
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2. Rotation 

For this case the system has enough energy to cover the full range of the angle q and 
there is no point where the momentum vanishes. The momentum is given by: 



p= \/2(E + a 2 ± cos q). (Ah 

The corresponding action integral is: 



i r^w 



27T./0 



J r = — / \ 2{E + af cos q) dq. (A9) 



Again using a change of variables we can transform this integral to: 

2«i r 



J r = — - / \/l-k 2 r sin 2 d<f>, (A10) 

kr^ JO 



where the modulus k r is given by k r = J "^ 2 ■ This integral is just twice the complete 
elliptic integral of the second kind, so that the action for the rotation case is: 

j r = ^l gl (k r ). (All) 

3. The separatrix case 

For the separatrix case we can repeat the same kind of calculations as for the rotation 
case. Noting that for this case we have k r — 1 and £i(l) = 1 the action for the separatrix 
case is just: 

J s = — . (A12) 

7T 

Appendix B: Action variables for the symmetric double- well 

In this appendix we compute action variables for the integrable symmetric double-well 
system. The Hamiltonian of the system is: 

In order to determine the action variable of this system we have to consider three different 
cases distinguished by the value of the energy E: E < 0, the separatrix case where E — 0, 
and E > 0. 
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1. Case E < 0. 

For this case the energy is below the saddle point energy, E = 0. Trajectories are confined 
to one of the two wells. Because the double- well is symmetric the action variable is the same 



for the motions in either well.From eq. (Bl ) we can write the momentum p as a function of 
the energy and q: 



,, v 2£ + a 2 g 2 -|. (B2) 



The associated action integral is: 



J 1 = ±-jpdq=±-2 J ^2E + a 2 q 2 -^dq, (B3) 

where a and b are the two roots of the equation p = for the well situated situated on the 
side where q > for example and with a < b. 

Using a change of variable we can rearrange the result into the form: 

J^^T f 9c "* 9 dO, (B4) 

v^Wo Vl " K sin 2 6 

where 5 = b 2 — a 2 and k\ = w|. The latter integral can be evaluated by parts and we 
finally obtain: 

h = ^^ [(2 - kj^ih) + (2kj - 2)/^)] . (B5) 

2. Case E > 0. 

For this case the energy is above the saddle point energy. In this case the equation 
p = has only two roots situated symmetrically with respect to the coordinate origin q = 0. 
Denoting the two roots a and b, a < b, the action integral is: 

1 I 4 f b 

J 3 = — (t> p dq = — / pdq, (B6) 



2tt J 2tt j,, 



Using an appropriate change of variable we can transform this integral to the form: 

25b 2 r /2 



Ti 







J 3 = / Jl-kl sin 2 9 sin 2 9 dO, (B7) 



where 5 = b 2 — a 2 and k% = wi. Inetgrals of this type are tabulated (see for example ref. 
I50|) and one gets the result: 



J 3 = ^Sr [(2*2 - l)£i(*a) + (1 - fc 3 2 )^i(*3)] • (B8) 



37rA;3 
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3. Case E = 0. 



The action in this case is obtained by setting a = 0, b = a/2«2, 5 = b 2 and k 3 = 1 in the 



previous case, and noting also that instead of having a factor 4 in eq. (B6) we have here a 



factor 2. Taking into account the fact that £i(l) = 1 we get for the separatrix action: 

3/2 



J * = it (B9) 
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Figure captions 



FIG. 1: 

product 
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product 
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OL\ = 0. 
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Phase portrait for example 1 (cf. sec. Ill A 1 ). The phase space is the Cartesian 



(denoted by "X" ) of the phase spaces for the uncoupled subsystems. 



Phase portrait for example 2 (cf. sec. Ill A 2). The phase space is the Cartesian 
(denoted by "X" ) of the phase spaces for the uncoupled subsystems. 



Dividing surface for example 1 (eq. (3.9)) for three different energies. Parameter 
8. (a) E = 0.5. (b) E = aj. (c) E = 1.0. 
. 4: Flux 4>(E) across the DS for example 1 as a function of energy E (a\ = 0.8). 



Dividing surface for example 2 (eq. (3.9)) for three different energies. Parameters 
= 1. (a) E = -0.15. (b) E = 0. (c) E = 1.0. 

Flux 4>(E) across the DS for example 2 as a function of energy E (a.\ — ol-i — 1). 
Gap time distributions, (a) Example 1, energy E? = 5. (b) Example 2, energy 
Note the existence of a singularity in the gap time distributions for both cases. 
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